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We perform a one-loop calculation of the strange quark polarization (As) of the nucleon in a SU(3) 
chiral potential model. We find that if the intermediate excited quark states are summed over in a 
proper way, i.e., summed up to a given energy instead of given radial and orbital quantum numbers, 
As turns out to be almost independent of all the model parameters: quark masses, scalar- and vector- 
potential strengths. The contribution from the quark-antiquark pair creation and annihilation ll Z" 
diagrams is found to be significant. Our numerical results agree quite reasonably with experiments 
and lattice QCD calculations. 
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The intrinsic strangeness content of the nucleon is a 
key ingredient to understand the structure and dynam- 
ics inside baryons. While the experimental investigation 
of the nucleon spin structure |lj clearly indicates that a 
strange quark sea exists and is also polarized relative to 
the nucleon spin, the successes of the naive spin-flavor 
SU(6) valence quark model in various aspects suggest 
that the strangeness content should belong to higher or- 
der effects for the nucleon. The SU(3) flavor chiral quark 
model, which couples light quarks to octet pseudoscalar 
mesons by the requirement of chiral symmetry, provides a 
natural mechanism for such a perturbative picture: at ze- 
roth order the ground state octet baryons are described 
by a SU(6) wave function of three valence quarks, and 
at second order in the quark-meson coupling sea quarks 
can be generated by emitting a meson from the valence 
quark. For example, in the nucleon the strange quark 
can be generated by emitting a K + from the u quark or 
by emitting a K° from the d quark (Fig. 1), and hence 
strange quarks can contribute to the nucleon spin struc- 
ture. 
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FIG. 1. Strange quark contribution to the nucleon spin; a 
cross on the quark line denotes the axial vector vertex 7 3 7 5 . 



In this paper we will adopt the standard perturba- 
tion theory to calculate the strange quark polarization 
of the nucleon in the framework of a SU(3) chiral po- 
tential model. As will be shown, up to second order the 
diagrams of Fig. 1 are the only contributions. Hence, a 
quantitative determination of the strange quark polariza- 
tion is a clean test for the interaction picture of the chiral 



quark model, which allows in this case only admixtures 
of K mesons. On the other hand the u or d quark polar- 
izations obtain contributions from many other diagrams 
and therefore depend on many more model parameters. 

To set up our calculation scheme, we first define the 
effective Lagrangian for the SU(3) chiral potential model: 

C = ip[ip- S(r) -j°V(r)]ip - 

-L^[5(r)(«r + * 7 5 AV0 + (a + z 7 5 A^)S(r)]V> + 
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The model Lagrangian is derived from the a model in 
which meson fields are introduced to restore chiral sym- 
metry jjj. ip is the quark field with flavor and color 
indices implied; the scalar term S(r) — cr + m repre- 
sents the the linear scalar confinement potential cr and 
the quark mass matrix m; V(r) — —a/r is the Coulomb 
type vector potential and F w =93MeV is the pion decay 
constant, a and 4>i (i runs from 1 to 8) are the scalar 
and pseudoscalar meson fields, respectively and \ are 
the Gell-Mann matrices. The quark-meson interaction 
term of Eq.(Q) is symmetrized since the mass matrix m 
does not commute with all for different quark masses. 
The zeroth order quark Hamiltonian is set up as 



H q = [ d 3 x^[a ■ -d + pS(r) + V(r)\ij>. 



(2) 



It has discrete eigenstates which are obtained by numeri- 
cal solution of the Dirac equation with a scalar and vector 
field We write the solution as: 
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Eq. (|3|) forms the basis of our unperturbed wave func- 
tions, where quarks are bound permanently by the con- 
finement potential which is included in H q . From Eq. (|^) 
we can construct the quark propagator: 
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D(x x ,x 2 ) = (0|T{^(n),^(a; 2 )}|0> 

= 9{t\ - t 2 ) ^ u a(xi)u a (x 2 ) ~ 

a 

O(t 2 -ti)^2v (xi)vj3(x 2 ), (4) 



The meson propagator given by Eq. (|l|) is the free one: 
&ij(xx,X2) = (0|T{^(a;i),^(a;2)}|0> 

A T _ 2 , ... ■ (5) 



(27T) 



q z — mf + it 



Given the unperturbed basis we can construct any 
physical quantity up to a desired order in the quark- 
meson interaction. In the following we are studying the 
quark contribution of flavor q (q = u, d, s) to the nucleon 
spin which is defined through 
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At zeroth order H q gives the usual SU(6) three-quark 
states for the nucleon with the single quark wave func- 
tion u a in the ground state. The zeroth order diagram 
for the numerator of Eq. (Q) is indicated in Fig. 2A, and 
the denominator by the diagram of Fig. 3A, which is 
simply unity. Clearly, strange quarks do not contribute 
at this order. 

The corresponding Feynman diagrams which con- 
tribute to As up to second order are shown in Figs. 
2 and 3. The denominator (iV|iV) can be denoted as 
(1 + const. /F%), which can be expanded (1 — const./ F% + 
• • •) and has then to be multiplied with the numerator 
(N\ J d 3 £i/>7 3 7 5 ?/;|iV}. If finally, only terms of order 1/F% 
are kept in the product of the normalization and the ma- 
trix element of the spin, this has no effect on As, since 
already the lowest order admixture of s quark is propor- 
tional to 1/F%. 
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FIG. 2. Feynman diagrams for the matrix element 
(iV| J d 3 xtp'y 3 'y 5 tp\N} up to second order; A is of the zeroth 
order, B is the renormalization counter term, C and D are 
vertex and exchange diagrams respectively. The meson line 
in C can be a ir and f] (while the intermediate quark is it or 
d), or a K (while the intermediate quark is s); the meson 
line in D can only be a -n or r\. 
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FIG. 3. Feynman diagrams for the normalization (N\N) 
up to second order; A is of the zeroth order which is simply 
unity, B is the meson exchange diagram. The meson line in 
B is a 7r or r\. 



In the Lagrangian of Eq. (|JJ) the main effect of the 
nonperturbative quark-gluon interaction is supposed to 
be included by the scalar and vector potentials. In prin- 
ciple we can also include a residual perturbative gluon 
piece. This will introduce further modifications on Au 
and Ad. However since the perturbative quark-gluon in- 
teraction is diagonal in flavor space, it cannot generate 
strange quark admixtures in second order for the nucleon. 

Now we are in a good position to calculate As for 
the nucleon: up to second order the only diagrams we 
need to consider are the subset of the diagrams of Fig. 
2, which are given in Fig. 1. For the evaluation, we 
first give the explicit form for u a {x) and vp{x) with: 
u a {x) = e~ lEai u a {x)r a , vp(x) = e %Eflt vp{x)T^, where 
r is the flavor wavefunction and the spatial wavefunction 
is: 



u a {x) 



9njl 

'i^ ■ rfnji 
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where g and / are real functions, n is the radial quantum 
number, and F^™(r^) are the vector spherical harmonics. 
The second equality of Eq.(j7j) follows from a ■ rYJp = 
— Yjjt with /' = 2j — I. For computational convenience, 
we will use exactly the same form for v$ (x) . Since for the 
antiquark solution the lower component is the large com- 
ponent, for vp(x) I is actually the orbital quantum num- 
ber of the small component, and |-Ej=j+i/a| > \Ej = i_ \/ 2 \. 
Thus for the antiquarks the sequence is inversed. 

Denoting the initial and final quark states as Ui and 
u f respectively, the contribution of the diagrams of Fig. 
I is: 



5s = i / d :i xd' i x 1 d' i x 2 A(x 2 ,xi)uf{x 2 )S{r 2 )^\ 1 x 



D(x 2 ,x)YrD(x, x^Sinh^X'uiixx) 



(8) 



Here we use Ss to indicate that it is only the contribution 
from a single quark state. Inserting the explicit expres- 
sions for the propagators, we get 
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= -j^ I d 4 x 1 d 4 x 2 u f (x 2 )S(r 2 )~f 5 \ l x 



6(t 2 - t)6(t - t\) ^2 u a {x 2 )A aal u a >(xi)- 
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9(t! - t)9(t - t 2 )J2vp(x 2 )App / v p ,(x 1 ) - 

W 

9(t 2 -f)0(*i - t) };U a (x 2 )A al3 ,V0, (xi) 

a/3' 

Vp(x 2 )Ai3 a ,U a >(Xi) 

/3a' 

SinrfXu^)-— / 5—^. (9) 
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q 2 — rn~ + ie 



where A aa > = J d 3 xu a j 3 ^ 5 u a ' , and similarly for A^ 
etc. The four time-ordered terms in Eq.(^|) correspond 
to the time-ordered diagrams of Fig. 4. 
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FIG. 4. Time-ordered diagrams of Fig. 1; A is the posi- 
tive-energy state contribution; B is the negative-energy state 
contribution; C and D are the quark-antiquark pair creation 
and annihilation U Z" diagrams. 

We omit here the details for calculating 6s of Eq. <^). 
The integrals of Eq. (||) can be reduced analytically to 
radial integrations at the vertex points (r\ and r 2 ) and 
of loop momentum The the remaining integrations 
are carried out numerically. As for the whole nucleon is 
just 6s times a spin-isospin factor which can be straight- 
forwardly calculated to be 2. 

In Table I we list our model parameters. Since = 
93MeV and vtlk — 495Mev are fixed by experiment, 
our model contains four free parameters: the two quark 
masses m u ^, m s and the two strength constants of the 
scalar and vector potential denoted by c and a. The 
parameter a is fixed by the long-wavelength, transverse 
fluctuations of the QCD based static-source flux-tube pic- 
ture It was obtained to be 0.26 in || and 0.30 in 
M, while a much larger value of about 0.52 was used 
by the Cornell group Recent lattice calculation || 
got a value around 0.32 in the quenched approximation, 
and suggested that relaxing the quenched approximation 
may lead to a ~ 0.40. Quark masses and confinement 
strength are rather uncertain quantities. To study the 
variation of As over all the parameters, we choose in our 
calculation four different sets of parameters, including 
both current and constituent quark masses. 

We study very different parameters because they often 
vary significantly from one model to another. For this 



model one possible choice of parameters to produce the 
correct nucleon mass, <?a, etc. is given in Ref. |Tc|| . 

Table II gives the numerical results of As for 
the first two sets of parameters. The intermediate 
quark/antiquark states are summed over up to a radial 
quantum number of n = 8 and total angular momentum 
j = 17/2. We also list the intermediate results with the 
summation including states up to n = 6 and j = 11/2. 
The contributions from the four time-ordered diagrams 
in Fig. 4 are given separately. We note significant con- 
tributions from Fig. 4C and Fig. 4D, in which a quark- 
antiquark pair is created or annihilated by the axial vec- 
tor current; these processes are usually referred to as the 
"Z" diagrams. On the other hand the diagram of Fig. 
4B gives a fairly large positive contribution, therefore 
if the "Z" diagrams are neglected we would incorrectly 
conclude that As in the nucleon is positive. 

From Table II one would conclude that a stronger con- 
finement also gives a larger As. This is due to the cou- 
pling of the meson field to the quark field which is pro- 
portional to the effective quark mass S(r) — cr + m. 
However, to compare with the energy scale in the lat- 
tice QCD calculation of As, we should sum the excited 
states up to a given energy instead of given radial and 
orbital quantum number. The resummed As according 
to energy are given in Fig. 5. Since the quark states are 
discrete, we get plateaus in Fig. 5 at the energies where 
no new states emerge. 

TABLE I. Model parameters 



para. 


m u ,d 


m s 


Q 


c 


set 


[MeV] 


[MeV] 




[GeV 2 ] 


1 


10 


150 


0.26 


0.11 


2 


10 


150 


0.26 


0.16 


3 


300 


500 


0.26 


0.11 


4 


10 


150 


0.50 


0.18 



TABLE II. Numerical results for As by summing the inter- 
mediate quark/antiquark states up to given radial and orbital 
quantum numbers. 
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set Fig. 4A Fig. 4B Fig. 4C Fig. 4D sum 



6 11/2 



.0212 
.0340 



+.0829 
+.1238 



.0844 
.1255 



.0844 -.1070 
.1255 -.1613 



17/2 1 -.0220 +.0964 -.1109 -.1109 -.1475 
2 -.0345 +.1445 -.1655 -.1655 -.2211 



TABLE III. Numerical results for As by summing over the 
intermediate quark/antiquark states up to the energy 1.7GeV. 



para, set 


Fig. 4A 


Fig. 4B 


Fig. 4C 


Fig. 4D 


sum 


1 


-.0228 


+.0785 


-.0742 


-.0742 


-.0927 


2 


-.0365 


+.0958 


-.0805 


-.0805 


-.1017 


3 


-.0464 


+.1538 


-.1046 


-.1046 


-.1018 


4 


-.0400 


+.0834 


-.0636 


-.0636 


-.0838 
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Since the strange axial current is a non-conserved com- 
posite operator, it has divergent matrix element (as is 
seen in Fig. 5), and therefore must be renormalized. 
Analogous to the lattice renormalization, we cut the 
quark intermediate states at an energy of 1.7GeV, which 
is roughly the inverse of the lattice spacing in the lattice 
calculation of As (a -1 = 1.74GeV in The "renor- 

malized" results are given in Tabic III. (In principle, we 
can also do renormalization by imposing a cutoff on the 
meson momentum, such as using the Pauli-Villars regu- 
lator (q 2 — A 2 ) -1 . But then the quark intermediate states 
have to be summed up to convergence. In practice this is 
not workable. An illustration with A = 1.7GeV is given 
in Fig. 5.) 
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FIG. 5. plot of As as a function of the maximal energy 
up to which the intermediate states are summed. 

We note a very interesting phenomena in Fig. 5: the 
result for As summed up to a given energy is rather ro- 
bust against the variation of all the parameters. The in- 
sensitivity is especially impressive compared to the huge 
variation of m and c. Table III shows that the "Z" dia- 
gram's contribution is still significant. 

The insensitivity of As on the parameter sets can be 
attributed to the fact that the increase of m and c (see 
Eq. ([!])) enhances the quark-meson coupling and moves 
up the single quark state energy. Thus the contribution 
from a single quark state increases due to the stronger 
coupling but less states are accessible to be summed over 
up to a given energy Similarly, the increase of a sup- 
presses the contribution of a single quark state since the 
lower Dirac components of the quark wavefunctions are 
increasing. But it also reduces the quark state energy, so 
we have more states to sum over. 

The main results of this paper can be summarized as 
follows: 1) Strange quark polarization is a very clean and 
robust prediction of the chiral potential model. Up to sec- 
ond order the only contribution arises from the diagram 
of Fig. 1. As depends only weakly on the model parame- 
ters, and our calculation shows further that the variation 



of these parameters does not influence As too much, pro- 
vided we sum over the intermediate quark state up to a 
given energy. 2) The contribution from the intermediate 
excited quark states are important. It is not enough to 
restrict the intermediate state to the ground or the first 
few states. 3) Among the time-ordered diagrams, the 
quark-antiquark pair creation and annihilation "Z" dia- 
grams are significant. It is the "Z" diagrams (Figs. 4C 
and 4D) that introduce a negative value for As in the 
nucleon, while the intermediate negative-energy states 
(Fig. 4B) gives a fairly large positive contribution. The 
importance of the pair creation and annihilation contri- 
bution to Aq has also been noticed by some of us previ- 
ously in a valence and sea quark mixing model [ jyj . 4) 
Our numerical result is quite consistent with experiments 
(As(Q 2 = 3GeV 2 )= -0.10±0.01±.- ■, where the second 
i sign represents further sources of error, principally the 
low x extrapolation |l^]) and lattice QCD calculations 
(As = -0.12(1) |§, -0.109(30) @), and is also consis- 
tent with a schematic calculation in the context of chiral 
quark model by Cheng and Li |1| . To the best of our 
knowledge this is first time that As is consistently calcu- 
lated up to the one-loop level in a quark model. 
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